Effect of inter-particle interaction in a free oscillation atomic interferometer 



Thomas Fogarty 1,2 , Anthony Kiely 2 , Steve Campbell 1,2 and Thomas Busch 1,2 
1 Quantum Systems Unit, OlST Graduate University, Okinawa 904-0495, Japan 
2 Department of Physics, University College Cork, Republic of Ireland 
(Dated: April 30, 2013) 

We investigate the dynamics of two interacting bosons repeatedly scattering off a beam-splitter in a free 
oscillation atom interferometer. Using the inter-particle scattering length and the beam splitter probabilites 
as our control parameters, we show that even in a simple setup like this a wide range of strongly correlated 
quantum states can be created. This in particular includes the NOON state, which maximizes the quantum 
Fisher information and is a foremost state in quantum metrology. 
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I. INTRODUCTION 

In all scientific pursuits accurate measurements are cru- 
cial and many of the most successful techniques in quantum 
metrology make use of the principles of interferometry. Re- 
cently significant progress has been made by recognizing that 
the use of quantum correlations, in particular entanglement, 
enables us to make the most precise measurements physically 
allowed by the Heisenberg uncertainty principle d-ID. Opti- 
cal interferometers are able to generate a wide range of quan- 
tum correlated states, such as the NOON state [5 - 8 ], however 
a major draw back are their short coherence times. In order to 
enhance measurement precision through the use of entangle- 
ment, longer coherence times are highly desirable. 

Making use of atomic ensembles can enhance the lifetime 
of a generated state, however the often unwanted and hard- 
to-control scattering interactions can make them difficult to 
work with [2]. Nevertheless, some remarkable progress has 
been made using Bose-Einstein condensates (BECs) as re- 
sources [9l 02). For example, in the presence of attrac- 
tive interactions these allow for the formation of bright soli- 
ton states, which are non-dispersive and have been suggested 
as good candidates for the creation of macroscopic spatial su- 
perpositions fTTlfT4ll . Furthermore, the ubiquitous presence 
of harmonic traps for ultracold atoms has led to new ideas 
for interferometry designs based on the periodic trap dynam- 
ics lfT5HT8l . Such schemes, which present a viable approach 
to atomic interferometry require often minimal experimental 
efforts and are referred to as free oscillation atom interferom- 
eters. 

Here we investigate the behaviour of two ultracold atoms in 
such a free oscillation atom interferometer and fully take their 
mutual interaction into account. We start with two bosonic 
atoms located on one side of a harmonic trap split by a delta- 
potential, whose strength can be adjusted. The atoms are then 
released and allowed to scatter off the barrier twice, thus real- 
izing a Michelson type interferometer. By employing numer- 
ical diagonalization techniques we are able to exactly solve 
the model and determine the atom pair's full density matrix 
at any moment in time. While previous studies have explored 
how different properties of the trap affect the performance of 
an interferometer lfl~5l [T9l , here we rigorously assess the ef- 
fects which different interaction regimes and beam splitter ra- 
tios have on the non-classical nature of the states created. We 



quantify this by calculating the Quantum Fisher information 
(QFI) [20], and show that for a certain range of parameters this 
simple setup can generate the highly desirable NOON state. 

The remainder of the paper is organized as follows. In 
Sec. [Il] we formalize the physical model and present the var- 
ious tools to be used throughout. In Sec. Ill A we assess the 
case when the atoms possess an attractive interaction, while in 
Sec. |IIIB| we explore the repulsive regime. The experimental 
feasibility is considered in Sec. [IV] and in Sec. [V] we present 
our conclusions and discussions of the results. 



II. PRELIMINARIES 



A. The Model 

The atomic interferometer we consider is a harmonic trap 
punctuated centrally by a delta function potential. The delta 
function barrier will act as a beam splitter for the interacting 
atoms, and for simplicity we restrict our investigation to the 
case of two atoms. We assume the trap is such that only lon- 
gitudinal motion is permitted and transverse motion is tightly 
restricted, thus forming a quasi one-dimensional system. The 
Hamiltonian is then given by 

n- 1 \ " / 

(1) 

where m is the mass of each particle, Q the frequency of the 
harmonic potential and kq is the height of the ^-function bar- 
rier. Throughout the paper, unless otherwise stated, all units 
are dimensionless. At low temperatures the boson-boson in- 
teraction, V, can be approximated by a point-like potential 



V (|*i - x 2 \) = gw S Qxi -x 2 \), 



(2) 



where g\ D is the one-dimensional coupling constant between 
particles defined in terms of the three-dimensional scattering 



length as g 1D 
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with C 1.4603 and a ± 



112111 . This parameter will be central in our analysis of different 
regimes and can be experimentally tuned by applying a Fesh- 
bach resonance, a powerful technique that is well established 
in cold atomic physics ll22l . 
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and must be solved numerically using, for example, a discrete 
variable representation (DVR) method El |23. The DVR 
method allows exact diagonalization of the many body Hamil- 
tonian and scales as N P N , where N p is the number of points 
taken in configuration space. In general this is numerically 
expensive, however the restriction of our analysis to N = 2 
particles allows the calculations to be tractable. Time evolu- 
tion is then achieved by constructing the time dependent wave 
function in terms of the eigenstates of the hamiltonian < Kq 
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FIG. 1: (Color Online). Single particle density versus time for (a) 
two attractive and (b) two repulsive atoms. The barrier is positioned 
at x = and the two particles are initially trapped at d = 6 with 
e = 5.164. At time t s the particles scatter off the barrier and come 
to rest at time t A at the classical turning point. At time 3t s the atoms 
recombine and scatter a second time and come to rest again at time 

tB- 



Initially the two atoms are prepared in a separate tight har- 
monic trap a distance d from the centre of the interferometer 
trap and their state is given by the ground state of the Hamil- 
tonian 



in which 



2 

n=\ 



h 2 d 2 



1 



O A 2 + ^> m0J ( Xn 

2m oxt 2 



df \+g' w {\xi-x 2 \). (3) 



Here co is the trap frequency of the preparatory trap and 
the interaction is modified to reflect this different environ- 
ment. In the following we will make use of natural units 
such that the coordinates are rescaled with respect to the char- 
acteristic scales of the harmonic oscillator, x n = x n ja^ and 
E n = E n /(HQ). Here a ± = ^ti/m£l is the width of the trap 
ground state. Thus we have 
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is the overlap of the individual solutions of the two Hamil- 
tonians. Due to the atom's initial potential energy they will 
gain velocity, scatter at the barrier at time t s = n/2£l s (scat- 
tering A) and return to the classical turning points of the trap 
at tA = n/Qs (see the dynamics of the single particle den- 
sity in Fig. [T]). Here Q«j < £1 is an effective trap frequency 
adjusted to the presence of the delta function barrier. At 
time 3t s = 3n/2£ls the atoms scatter a second time (scat- 
tering B) and again return to the classical turning points at 
ts = 2^/Qj. This setup resembles an atomic Michelson in- 
terferometer. While the following analysis can easily be per- 
formed by describing the barrier with a well localised poten- 
tial of any shape, our choice of a delta-function is done to 
clearly isolate the interesting physical effects and does not 
constitute any loss of generality. A delta- function potential is 
a good approximation to a localised laser potential or an inter- 
action with an atomic impurity fixed at x = 0. In the first case 
the barrier height k can be experimentally tuned by chang- 
ing the laser intensity, whereas in the second case a Feshbach 
resonance can be employed. This, coupled with the capacity 
to alter the inter-particle interaction, means we have a highly 
adaptable system with which to create superposition states. 



^ = Z ( _ h txi + \ e " J ) 2 ) + 8 6 (l ^ " * 2|) ' (5) 

where e = cj/CI is the ratio of the preparatory trap frequency 
to the one of the interferometer trap, g - gwa_L = g 1Z) <2±/£ 2 , 
d - d/a ± and k - Koa ± . 

In order to solve the Hamiltonians, and < ?4 > , we 
must determine the single particle eigenstates and associ- 
ated energies. For the preparatory stage the time-independent 
Schrodinger equation 

9<aj^(xu x 2 ) = EJ/ n (x u x 2 ) (6) 

can be treated by taking advantage of the separability of 
the Hamiltonian into centre of mass and relative coordi- 
nates, for which the solutions are well known l23l . How- 
ever, the Schrodinger equation for the interferometer poten- 
tial < Ha0 w (fi, x 2 ) = E n <p(x\,x 2 ) does not allow such a luxury 



B. Quantum Fisher Information and von Neumann Entropy 

In the following we will thoroughly explore the features 
of the states at the times t& and ^. For this we will use the 
Quantum Fisher Information l26lL which defines the amount 
of information an observable can obtain about an unknown 
parameter. A state with maximal QFI will allow the most sen- 
sitive measurements, for example of relative phases. For a 
pure state \ij/((p)) the QFI is defined as 

T Q = 4 (W(tp)\ <A'(^)> - I W(<p)\ ^))| 2 ) , (9) 

where \iff'((f)) = d\ij/((p)) I dip. For separable states the max- 
imum QFI is equal to the number of particles (or quanta), 
N, used in the interferometer, which corresponds to the shot- 
noise limit. However one can go beyond this limit by using 
entangled states which can yield a maximum QFI of N 2 , the 
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Heisenberg limit fT9ll . One particularly important class of 
states that reach this limit are the so-called NOON states, 

Wnoon)= ^(|A0|0> + |0>|AT», (10) 

which in our scheme corresponds to both particles being si- 
multaneously on the left-hand side (LHS) and on the right- 
hand side (RHS) of the barrier. Thus we are we are looking 
at the spatial correlations of the two atoms 1271 . Of course, 
NOON states are not the only interesting non-classical states 
to study in interf erometry, however as we are examining N = 2 
they are the most prominent and our study is in line with the 
optical counterparts recent state of the art experiments l28l . 
For larger systems, i.e. iV > 2 we expect the scheme to be 
extremely versatile. 

The scheme generates a pure bi-partite entangled state, and 
as such we will also make use of the von Neumann entropy 
(vNE) to quantify the entanglement of the atoms. It is defined 
by the entropy of the reduced single particle density matrix, p, 
as 

S(p) = Tr[p log 2 p] = J] Mo&Ai, (11) 

i 

where the A t are the eigenvalues of this matrix and defined by 
f p(x\,X2)Xi(x2)dx2 - AiXi(x\). Due to the required symme- 
try of the wave-function for identical bosons one must be care- 
ful when dealing with the vNE as an entanglement measure in 
certain situations as discussed in [29]. In our setup, however, 
the dynamical scattering process and constant interaction be- 
tween the particles ensures that any finite von Neumann en- 
tropy signals genuine entanglement. Since the vNE measures 



the total entanglement, and therefore accounts for both inter- 
particle and spatial entanglement, it can be expected to show 
some qualitative differences to the QFI. 

III. ANALYSIS OF DIFFERENT INTERACTION REGIMES 
A. Attractive Interactions 

1. Scattering A 

We first examine the state of an attractive dimer after scat- 
tering once off the delta barrier. In Fig. [2] (a) we plot the QFI 
as a function of attractive interaction strength, g, and barrier 
height k. The thick black line signifies the classical shot noise 
limit at Tq = N = 2, which is attainable for separable states. 
Interestingly we find that even for a weakly interacting dimer 
we can exceed this bound. As we increase the attractive in- 
teractions between the atoms the QFI increases to its maximal 
obtainable value of N 2 = 4 for a barrier height of k « 1 . In 
panel (b) we see that the behavior of the vNE is qualitatively 
in agreement, although more complex. The small scale details 
are due to the inter-particle interaction leading to a constantly 
varying inter particle entanglement, which is not captured in 
the calculation of the QFI. Looking at the transmission coeffi- 
cient T, Fig.[2](c), we see that the maximum QFI is achieved 
for symmetric splitting (T = 0.5). To confirm the state gen- 
erated in this situation is the NOON state -^=(|20> + |02», 
we show the various population coefficients in Figs. [2] (d)- 
(f). One can immediately see that the region in which the QFI 
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FIG. 3: (Color Online). Contour plots for the (a) QFI, (b) vNE, (c) transmission coefficient (T), and population coefficients for states (d) |20), 
(e) |11), and (f) |02) as a function of attractive interaction strength g and barrier height k at time t B for the same initial state as in Fig. [2] 



is maximized corresponds to states for which the |11) com- 
ponent is suppressed and the |20) and |02) components are 
equally populated. This can intuitively be understood by real- 
ising that the relatively strong attractive interaction within the 
dimer makes it hard to split the pair of atoms into one on the 
left and one on the right. In fact, the situation is analogous 
to the one of bright, atomic solitons, where it has been shown 
that macroscopic superposition states can be created by mov- 
ing an atomic soliton through a barrier of finite width lfTTl[T4ll . 



2. Scattering B 

After the second scattering process the dynamics becomes 
more complex for the attractive dimer. Examining the QFI, 
Fig. [3] (a), we see that even for weakly attractive particles we 
can attain Tq ~ 4 and as we increase the interaction strength 
we find the QFI peak at two values of the barrier height, k. 
The behavior of the vNE, panel (b), is qualitatively similar 
and is also mirrored in the transmission coefficient, T, panel 
(c). Once again, for T = 0.5 we see a maximum QFI. The 
most striking feature is clearly the intricate series of maxima 
appearing in all panels. This is due to the phases accumulated 
by the atoms at the beam splitter and when traveling along its 
two arms. For the case of non- symmetric splitting the differ- 
ent interaction energies of the particles lead to a difference in 
phase, which in turn leads to the observed interference fringes. 
We see the same qualitative behavior in the various popula- 
tion coefficients shown in Fig. [3] (d) -(f), where the maximum 
QFI again corresponds to a suppression of the |11) state and 



an equal population of the other two states. Interestingly the 
value of k which resulted in a maximum QFI for scattering A 
results in a minimum QFI for scattering B for the same value 
of the interaction strength. 

To get a deeper insight into the process, we show in Fig. [4] 
the densities of the four lowest lying eigenstates, \x% (i = 
0, 1, 2, 3), of the reduced single particle density matrix at times 
tA and ts and for g = -7 and k = 0.4. Fig. [4] (a) shows that 
at tA (where we have Fq = 2.035) three of the orbitals are 
located on the RHS of the trap and one on the LHS. This is a 
result of the large attractive interaction, which does not allow 
the particles to become spatially split by the barrier and thus 
the transmission coefficient is either T = or T = 1. Even 
though each orbital occupies one side of the trap only, this is 



(b) 





FIG. 4: Densities of four lowest lying eigenstates of the reduced sin- 
gle particle density matrix for the attractive dimer at times (a) t A and 
(b) t B for k = 0.4 and g = -7. At t A the lowest, the second excited 
and the third excited orbitals occupy the RHS of the trap and the first 
excited orbital is localised on the LHS. At t B the ground state and 
the second excited orbitals are on the RHS, whereas the other two 
orbitals are on the LHS. 
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FIG. 5: (a) The orbital occupation numbers versus time for attrac- 
tively interacting particles k = 0.4 and g = -7. The darkest line cor- 
responding to lowest energy orbital, with each progressively lighter 
shade representing the next higher energy orbital, (b) The QFI (solid) 
and vNE (dashed) versus time. The maximum F Q reached is N 2 in- 
dicating that a NOON state is created. 



not a NOON state as can be seen by looking at the orbital oc- 
cupation numbers (see Fig|5](a)). We find that at this point 
the lowest orbital has still the largest occupation number and 
higher lying ones have lower occupation. Fig.[4](b) shows the 
situation at t B and we find two orbitals occupying each side of 
the trap. The orbital occupation probabilities for the ground 
and the first excited state are degenerate after the scattering 
event B, which proves the NOON nature of the state and ex- 
plains the resulting T Q = 3.9998 * N 2 (see Fig(5](b)). The 
fact that the occupations are not exactly degenerate close to t B 
is due to the inter-particle interaction, which is reflected in the 
dynamics of the vNE, dashed line Fig|5](b). It displays two 
pronounced dips exactly at t& and t B , indicating a prominent 
change in the internal structure (due to the re-focussing at the 
classical turning point). Making the choice oj = Q leads to 
perfect degeneracy for all times after scattering B (and before 
the next scattering event). Note that the step- wise behavior of 
the QFI is due to its sensitivity to only spatial correlations and 
it therefore only changes during the scattering process, while 
the constant interaction between the atoms gives rise to the 
varying vNE. 
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FIG. 6: (Color Online). Contour plots for the (a) QFI and (b) vNE 
at time t A after scattering once off the barrier and the (c) QFI and (d) 
vNE at time t B after scattering twice off the barrier, as a function of 
repulsive interaction strength g and barrier height k. The initial state 
is the same as in the attractive case discussed above. 
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B. Repulsive Interactions 

1. Scattering A 

We now turn our attention to the case of repulsive interac- 
tion between the atoms. This regime gives rise to behaviors 
that do not promote the generation of spatial entanglement 
easily as the repulsive nature prefers a situation in which one 
atom occupies each side of the trap. Fig. [6] (a) shows that 
at time tA this set-up cannot produce states that outperform 
the best classically attainable states (since Tq < 2) for the 
whole range of parameter space considered. The maximum 
Tq-2 occurs for a barrier height k = 1.51 regardless of the 
interaction strength g, reaching the classical limit for a trans- 
mission coefficient of T = 0.5 (not shown). The vNE (Fig. [6] 
(b)) is also maximized for T = 0.5, reaching 0.9 for g = 4 



FIG. 7: The densities of the two lowest atomic orbitals of the reduced 
single particle density matrix at times t A [top panel] and t B [bottom 
panel] with a repulsive interaction with g = 1 and k = 2.1. At t A 
the each orbital has equal probability to be in the LHS or RHS of the 
trap. At t B each orbital is at opposite sides of the trap, indicating a 
highly entangled quasi-NOON state. 



and increasing to a maximum of approximately 1 for strongly 
repulsive atoms. 



2. Scattering B 

Similar to the case of attractive interactions, the variety 
of states created after the second scattering process becomes 
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FIG. 8: (a) The orbital occupation numbers are plotted versus time 
for a repulsive interaction with g = 1 and k = 2.1. The darkest 
line corresponding to lowest energy orbital, with each progressively 
lighter shade representing the next higher energy orbital. At time t B 
the ground and first excited orbitals become nearly degenerate, while 
the higher lying ones are only weakly populated, (b) The behavior 
of the QFI (solid) and vNE (dashed) as a function of time. The max- 
imum F Q reached is 3.883. 



much richer due to the phase acquired by the various compo- 
nents of the two-particle state. States with a Tq > 2 can now 
be generated, however they are restricted to a much smaller 
area of the parameter space compared to the attractive inter- 
action. In Fig.[6](c) we see for that small repulsive interaction 
a QFI of Tq > 3.5 can be reached and the vNE in Fig.|6](d) 
shows qualitatively similar behavior. As the interaction g is 
increased, the atoms enter the Tonks-Girardeau (TG) regime 
and the QFI approaches its classical limit of 2, corresponding 
to the state \if/) = - |20> + -4 1 11) + \ |02> resulting from a 
50/50 splitting. 

The state with the maximum QFI in this regime is achieved 
for g = 1 and k = 2.1 and we show the corresponding lowest 
two eigenstates of the reduced single particle density matrix in 
Fig. [7] At tA (upper panel) each orbital occupies both sides of 
the trap with nearly equal probability due to the 50/50 split- 
ting of the barrier and corresponds to Fq « 2. At time ts 
(lower panel) each orbital almost fully localizes on one side 
of the trap and the respective occupation numbers approach 
double degeneracy, cf. Fig. [8] (a), indicating the presence of 
a superposition state in accordance with Tq - 3.883, Fig. [8] 
(b). Due to the relatively weak repulsive interaction strength 
we see the vNE grow monotonically. 



IV. EXPERIMENTAL REALIZATION 

Evidently the scheme presented here has an immediate ex- 
perimental appeal as many of its components are readily im- 
plementable. When a coherent bilocalised state is created the 
detection of this state can be achieved by measuring the fringe 
visibility of the two particle interference which is maximal 
in the presence of a NOON state lfl3l . This can be done by 
exploiting the free oscillations in the harmonic trap after re- 
moving the barrier and a simulation of these fringes is shown 
in Fig. [9] The solid line shows the pattern associated with 
the generated NOON state for F Q = 3.9998 at g = -7 and 
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FIG. 9: The beamsplitter is removed from the interferometer and the 
bi-localized system is allowed to recombine. Due to the coherent 
superposition an interference pattern is observed for F Q = 2.0023 
(red dashed line) and F Q = 3.9998 (black solid line). The difference 
in fringe contrast is apparent. 



k = 0.4 and the dashed line shows what one would obtain 
for a state near the shot noise limit, Fq = 2.0023 at g = -7 
and k = 1.0985. The difference in fringe contrast near the 
shot noise limit and near the Heisenberg limit can be clearly 
seen. One could also implement the scheme described in l27l 
where after the two sides of the trap are allowed to interfere at 
a beamsplitter the correlations are measured by counting the 
atoms collected at two different detectors. 



V. CONCLUSIONS AND DISCUSSIONS 

We have presented a comprehensive analysis of two inter- 
acting particles in a harmonic oscillator interferometer. By 
considering a wide range of parameters we have demonstrated 
the importance of the inter-particle interaction and its neces- 
sity in creating metrologically useful states. By employing ex- 
act numerical diagonalization methods we were able to study 
the type of states dynamically created and assess their value 
by studying the correlations via the Quantum Fisher Informa- 
tion and von Neumann entropy. The QFI is a useful met- 
ric for determining a states use in metrology and we found 
that the maximally achievable values depends strongly on the 
number of scattering events. After scattering on the barrier 
once, the attractively interacting particles were able to exceed 
the shot noise limit and even create NOON states for certain 
parameters. However, for repulsively interacting particles a 
single beam- splitting process does not allow to exceed this 
limit. After a second scattering from the barrier, thus realiz- 
ing a Michelson interferometer, we found that NOON states 
could be created for both kinds of interactions, even though 
the range of potential parameters in the repulsive case was 
more limited. As previously noted, although our study explic- 
itly considers a delta- function barrier, the same results hold if 
one replaces it with a Gaussian barrier of finite width. In this 
instance the exact values of interaction strength and barrier 
height for optimal state generation will be slightly different to 
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those found here, however the qualitative conclusions remain 
unaffected. 

Let us finally comment extending the presented results to 
larger particle numbers. Treating the two-particle system has 
allowed us to rigorously assess the effect the inter-particle in- 
teraction has on the generation of metrologically useful states, 
while also allowing us to explore the entanglement dynam- 
ics via the von Neumann entropy. Going beyond two parti- 
cles is computationally extremely costly, and it is clear from 
our study that for ensembles of repulsive atoms highly corre- 
lated states will be very sensitive to the parameters involved. 
In fact, first calculations in the TG regime have shown that 
the generation of NOON states cannot be achieved this way. 
However for attractive interactions, even for a larger number 
of particles in the system, the scheme presented here should 
realize NOON states. As we have shown, when the atoms are 
strongly attractive, the bonds formed between the atoms are 



extremely hard to break, hence after the scattering processes 
they are much more likely to remain spatially close. This is 
somewhat analogous to the behaviour of bright solitons re- 
cently analyzed 11511151 . 
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